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1. ISTROOCCTION 

The objective of the Grant NSG 3307 f froa the NASA Lewis Research 
Center to the School of Aerospace Engineering at Georgia Tech, was to 
develop a new numerical approach for computing unsteady compressible 
viscous flows. This approcch offers the capability of confining the region 
of coiq>utation to the viscous region of the flow. The viscous region is 
defined as the region where the vorticity is nonnegligible and the 
difference in dilatation between the potential flow and the real flow 
around the same geometry is also nonnegligible. The method was developed 
and tested. Also, an application of the procedure to the solution of the 
steady Navier-Stokes equations for incoaipressible internal flows is 
presented. 
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2. MATHEMATICAL FORMIILATION 


In this chapter the MthesMtlcal relationships that govern the tiae 
history of c:aipressible viscous flow around an arbitrary body are 
presented. The kineiatic and kinetic boundary conditions are stated and a 
■ethod for estinating the surface vorticity is described. In order to 
iaprove readability, the derivations of some of the equations used in this 
chapter are omitted in the suiin text and given in the appencices. 


Governing Equations 

In the absence of body forces, the Navier-Stokes equations for a 
coaq>ressible fluid with density P, viscosity v, thermal conductivity k and 
ratio of specific heat coefficients Y , may be vrritten in an inertial 
coordinate system as follows 


3V 

P — ♦ P(V.7)V 


“7p + 7.T 


( 2 . 1 ) 


The equation of continuity is given by 

II * 7.(PV) » 0 


( 2 . 2 ) 


The energy equation and the equation of state are given by 
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P - If ♦ p V.(^) - V.Vp • -7.q ♦ (2.3) 

p • PRT • (~)Ph (2.4) 

H«re^V:T reprcsenta Che dissipacion function given in cartesian coor** 
dinetea by 


**■ 1 

^:t ■ t . . 5^ 

ij 3x 


where is Che shear stress tensor and can be expressed as follows; 


T. . - 2m e. . ♦ . X(7.V) 

ij ij ij 


Here e^j is Che rate of strain tensor and X is the second coeffi* 


dent of viscosity. 


SCineaatics Expressed in Integrsl Representation 


The vorticicy 4s and the dilatation 6 are related to the velocity 
field V by: 


^ X V ■ 3 (2.5) 

7 . V • 6 (2.6) 


The kinematics of the problem, governed by equations (2.5) and (2.6) 
are elliptic in nature, requiring the specification of Neuman, Dirichlet or 
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■ixad typt of bou'idory conditions for velocity. These conditions ere 
required both st infinity, known ss the farstream condition, and on the 
solid surface. In the present study, Dirichlet type boundary conditions 
are prescribed on ^e boundaries. 

The velocity boundary conditions are 


V - V.. 


on b 


(2.7) 


and R is the fluid doaain bounded by the boundary b. For external flow 
probleas, the boundary b consists of the farstream boundary c and the body 
surface s. On s, the no^^slip condition is used, while V ■ is prescribed 
on c. 

•e> 

since the kinematic relationships are linear in V, the velocity 

vector V can be decomposed into a solenoidal part and an irrotational 
> 

part V 2 with homogeneous boundary condition for the latter. Thus 


with 



V . - 0 

V X Vj^ ■ u) 


having boundary conditions 



0 




( 2 . 8 ) 


(2.9a) 

(2.9b) 


(2.10a) 


(2.10b) 
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«nd, for V^i 


^ . Vj - S 

^ X Vj ■ 0 


with boundary conditions 



'2c 


- 0 


(2.1U) 

(2.11b) 


(2.12a) 

(2.12b) 


Ml ( 1 ) has shovn tl-at it is possible Co recast the kioenatic aspect 
of the ptoblea into an integral representation for the velocity V in tersw 
of the vorticity ui and the dilatation 6, For two dimensional flows this 
means 


V(r;,t) 



w?(f-f ) ♦8(f-f ) 

o o 



2 


dR ♦ V 

OB 


(2.13) 


In (2.13) R is the region where the vorticity and dilatation are non- 
negligible. At high and ooderate Reynolds number the dilatation is 
significant at distances from the body where Che vorticity is already 
negligible. Bence the approach does not seem as advantageous as in the 
incompressible case. Howeveri by using the potential flow solution around 
the same body, it will be shown in the next section that the domain of Che 
computations can be reduced to include only the region tdtere the vorticity 
and Che difference in dilatation between Che viscous flow and potential 
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flow are both non-negligible. 

Dm of th« PQfoti«l Flow Solution to Rtduct th» Dowin of Cowpuf tiono 
EqiMtion (2*13) can be wriccen in the following form 





aix(r-r ) ♦6(r-r ) 
o o 



dR 


1 

■ r J 


B(r-ro) . 

72 


R2 


(2.14) 


where R1 is the region of the flow where vorticity is non-negligible and R2 
is the rest of the domein* extending to infinity for external flow 
problems. This expression can be written for the potential flow around the 
same body as follows: 


p o 



B (r-r )dR 

P o 



2 



8 (r-r )dR 

-E 2__ 



s 


Y 

-E 


x(r-r )dS 
0 

!'*■ |2 



(2.15) 


where the subscript p indicates potential flow and Y is the vortex sheet 

P 

strength on the surface, s, of the body <hie to the potential flow. Because 
equation (2.15) is a limiting case of the general viscous compressible flow 
relation (2.14), it gives the potential velocity everywhere except at the 
surface where the equation is identically zero* 
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If rtgion R1 «xtandt far enough from the body, Che combioation of 
(2«14) and (2* IS) givee 




-fj 


wx(r-r ) 
0 


♦(6-B)(r-r ) 

g ■■ ^ dR 


RI 





s 


Y x(r-r )d* 

_£ 2 ___ 



(2.16) 


The relation (2.16) implies the following. First, (^~6p) in region 
R2 is aaall enough so that its effect on the velocity in region Rl is 
negligible (a detailed discussion of this aspect is given in Reference 11). 
Second, one needs to solve only in region Rl, which is a smaller region 
chan the region ^ere B is significant. 

The kinematic boundary condition for the external flow problem 
requires Chat Che velocity has to reach the freestream velocity at an 
infinite distance away from the solid surfaces. This requirement is 
referred to in this work as the farfield boundary conditions. This 
requirement is satisfied by equation (2.16). However, if a finice‘>dif~ 
ference method is used without any coordinate transformation, this bound* 
ary condition is difficult Co satisfy since the computational dosiain Co be 
included becomes very large. 

Kinetics - The Vorticity, Dilatation, Density, 
and Energy Transport Equations 

By taking Che curl of equation (2.1) and using equation (2.4) one 


obtains 
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- 7x V X w ♦ (^~) ^InP x5h ♦ Vx 7. tJ- (2.17) 


Siailxrly, taking the divargenca of aquation (2.1) and uaing aqua- 
tion (2.4) raaulta in 


II - -7. [(V.7)vj -(^)A -(~)h7^1ap -(^)^h.^nP 


♦ 7.(i V.?) 


(2.18) 


Spacialixing tha aquations for oro-diaanaional caaai raarranging 
tha taxM in aquations (2.17) and (2.18), and collacting tha coafficianta 
of 6 and u, aquationa (2.17) and (2.18) baeooM 


~ - 7.(Vu>) ♦ ♦ 4»(p,6,w,h) 


(2.19) 


II • - 7.(fe) ♦ ♦ X(P,e,w, h) 


( 2 . 20 ) 


Tha full dataila of tha darivation of the above equations are given 
in appendix B. 

Tha density and energy equations can be also written in terms of the 
derived variables (a and 6 as follows 


- -5l(vinP) - r 


(2.21) 




( 2 , 22 ) 
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«^«r« (I) 4tnoc«a ch« Mcoitudc of the vorcicicy vector. The tervs 9 t X* ^ 

tad • look like source cenu and ere given in appendix B. 

The governing equations (2.19) to (2.22) have been noo«disien* 

sionalised fay nonsalixing the variables with respect to the following 

reference quantities: distance, L; velocity, V ; density, o *, enthalpy, 

2 

V I and tine, L/V , where L is the characteristic length of the body. 

00 fl® 

This type of noraalisatioo leads to the following non~diaensional para* 
■eccrs: Mach nusiber, M ( Reynolds number. Re; and Prandtl number, Pr. 

€3 

Formulation for the Study of an jseHi sively Started Airfoil 


In this section, the mathematical forsailation discussed above is 
specialised and applied to the study of cos^ressible lasdnar flow past an 
impulsively started airfoil. The airfoil geometry and the grid system are 
generated through a conformal transformation which t insforms the airfoil 
into a unit circle. The airfoil chosen for the numerical study is a 
modified 9X Joukowski airfoil. 

By using the transformation relations given in appendix , the 
governing equations (2.19*2.22) are vnritten in the transformed plane and in 
a conservation form as shown below. 

The vorticity transport equation is 


3u) 

at 


b • b ‘sc’ '’'-'b* 


(2.23) 


where H is the scale factor. 

The dilatation transport equation is 


10 




( 2 . 24 ) 


Iht dtnsity trMtporc «qu«cion is 

r • j (-»•<«“»') ■ ^ 


(2. 23) 


1h« snsrgjf transport aquation is 


II . K (-».(«,)) . ij » e 


( 2 . 26 ) 


Ohara ^ is tha apparant valocity in tha transfomad plana and tha divar* 
ganca and Laplacian oparators ara applied in tha transforaad plana. Tha 
sourca tarns 4tX 'I* and 6 are given in appendix ‘A. 

Tha radial and tangential conponents of tha valocity in the trans- 
fomad plane ara 


V' 




-i- ( 

•Iei; 

2uJ 


0 


-i- 1 

r <oH 

2ii J 


R1 

_ 1 


2iJ- 


Y .H(rcos(9-6 )-r )ds 


o 0 




♦ V 


e 


o o 




■h! 


H (3-8 )r‘«in(«-« )dR 


R1 


i‘« ♦ 1 2 


( 2 . 27 ) 


Ir 

la 0 1 


k/ 


%i 


u)H r sin(0-6 )dR 
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R1 


(S-6 Kr cot(e-0 )-r )dR 

P 0 0 


♦ V 


(2.28) 


Bquatioat (2.27) and (2.28) are essentially the same as the two 
conpoaents of equation (2.16) in cylinderical coordinates, except for the 
scale factor H. 


Surface Vorticity Deteraination 

The vorticity values avay fron the surface are deterained using the 
vorticity transport equation (2.23). In order to solve this equation, it 
is necessary to prescribe the vorticity values on the solid surface at all 
tine levels. To do that, the viscous region is conveniently divided into a 
vortex sheet of strength Y located on the surface, and an outer vorticity 
field «ihere the vorticity (ti and the dilatation 6 are assumed to be knotm. 
Applying equation (2.27) on the surface of the body yields 


- 0 


, /•H(Y“Y^)(r,cos(e-9 )-r^)d6 

I I p > o s 

8 ■*‘^8 "^*^8 
8 S 8 S O 


♦ V. 


(2.29) 


Here Y tepresents the integrated value of u)dn on the firet cell adjacent to 
the surface and r is the position vector for the points on the surface 
where the tangential component is calculated. is the tangential 

velocity at the body surface due to both the outer vorticity field and the 
whole dilatation field in Rl, and is given by 
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wH^(r co8(0-0 )-r )rdrd0 
o o 

2 2 

r r -2rr cos (0-0 ) 
o o o 



(B-8 )H^r 8in(0-0 )rdrd0 
P o 

-2rr coa(0-0 ) 
o o o 

(2.30) 


It oust be noted that the region R1 in the first integral does not 

include points on die solid surface. 

If r approach r , the first integral on the right-hand side of 
s 

equation (2.29) becones (20 j. 


/ (Y-Y.)H (r. cos(0-0^)-r )rd0 

P * 2 ! 

rj -2rf cos(0-0^) 

8 S S 5 O 


41T 


2ir 

1 ' 


(Y-Y„)Hds - i (y-Y„)H 
P 2 p 

(2.31) 


The principle of conservation of total vorticity gives 


o (Rl-s) 

and, since che solution is started by a non-*circulatory potential flow, it 
follows that 


27T 


4rr J 


Yp.H r^ d0 


(2.33) 


Substituting equations (2.30-2.33) into equation (2.29) yields 


1 f - 

J XT’ 


2 2 2 

-r^)(uH^ rdrde 


I f 


sin(0-0 )rdrd0 
o 


(Rl-s) '2rr^ cos(e-6^) cosv0-0^) 


♦ Y H 


(2.34) 
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In equation (2.34) Che region <Rl>s) is the compucacion region R1 
excludiiy Che surface a. Since Che radius of Che circle in Che Cransforacd 
plana is taken to be unity* equation (2.34) is rewritten as 


Y 


f 

ca8^(l-r^)rdrd0 

-- 1 

2im J 

l+r^-2r cos(0-0 ) 

ItH J 


(Rl-a) R1 


H^(S“6 )r sin(0-8 ;rdrdO 
P o 

l+r^-2r cos (0-0 ) 
o 



(2.35) 


Seg«entation of the Velocity Field 

The advantage of using equation (2.27) to calculate exterior flow 
problesK steas froa its ei^licit nature. Thus* the integral foraulation 
peraics the deteraination of the velocity on the boundaries of rectangular 
regions without regard to the interior nodes. In several cases, since 
rapid finite-difference computational schemes are available for solving 
the Poisson's equation in regions with rectangular boundaries, a combi- 
nation of equation (2.27) and such sdiemes can provide a faster way to 
compute velocities in exterior flow problems. For this reason, the 
computational domain is divided into compartments in which Che kinematic 
confutations are perfc^ed independently of each other. The choice of the 
scheme* to be applied in each compartment, dep«ids upon the shape of Che 
body surface and on Che relative distance between the body surface and the 
compartment. For exanple, as will be shown later in the static stall case, 
Che integral relation (2.27) is used in the whole wake and in inner regions 
adjacent Co the airfoil surface in order to compute the velocities. The 
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Poiason's equation is used in the rest of the computet ions 1 domain and the 
velocity on the boundaries is calculated using the integral relation (2.27). 

The Poisson's equation for the tangential velocity in the trans**’ 
formed ;>lane is derived as follows: 

The vorticity and dilatation can be written as 




X V ■ oil' 


!^..2 _ ■*: 


Ui 


(2.36) 


V * y • SR • B 
o 


(2.37) 


Upon taking the oirl of the terms in equation (2.36)t the equation becomes 


VxW »7x7xV» 7(7. V) - T^V 
o 


(2.38) 


Substituting equation (2.37) in equation (2.38) one obtains: 


T^V 


70 - 7 X w 
o o 


(2.39) 


Writing equation (2.39) in polar coordinates 


(l»“o j=“o 1 2 

*. - ■KS ■= — } e -(r — + - -r=- > e- * V 


3w- 


38 . 


7 V. e„ *7 V e = 0 (2.4C) 


/r 38 3r pr /3r ^ r "0 " ''0 "r «r 


By singling out the tangential component terms in equation (2.40) one gets 





1^ 
r 30 



(2.41) 


Equation (2.41) is the Poisson's equation for the tangential velo** 
city written in the transformed plane. 


Once the tangential component of the velocity, is determined the 
radial component of the velocity is calculated explicitly hy using the 
definition of the dilatation in die transformed plane, namely 


7 . V ■ Bh « 6 




6 

o 


(2.42) 


Pressure and Shear Calculations 

Since the surface pressure details are needed for any load esti- 
mation cn the body surface, the equation of state (2.4) is used to de- 
termine die pressure on the surface as follows: 


p * (-^— )Ph 


A pressure coefficient can be written as 




(2.43) 


(2.44) 


Mherc is a reference pressure. 

As will be shown later, the gradients of the flow variables on the 
upper surface of a stalled airfoil are very sensitive to small disturbances 
created either using different approximations to die governing 
equations or tqr adopting different boundary conditions. Compressibility 
effects are expected to be small for - 0.4, ^tnd in order 
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to copture these snail effects, it was decided Lo compare Che present 
compressible results with incompressible results obtained by using exactly 
Che same mathematical and numerical procedures. For this reason, a 
different (from that employed in the test cases) scheme for computing the 
surface pressure, similar to Che scheme used in the incompressible case 
(7 ), is developed and presented below. 

.In the body-fitted coordinate system, the vector momentum equation 

is 


^ at * * - ^p + 7. T 


(2.45) 


AC the surface, the acmentum equation is reduced to the following 
timple form because of Che no-slip condition. 


7P » 7. T 


(2.46) 


Talcing the dot product of the above equation with the tangential 
unit vector t at the surface, defined positive in the counterclockwise 
sense one gets 




3s 


(7.T).t 


(2.47) 


where s is the coordinate direction tangential to the surface, and is 


measured positive in the counterclockwise sense. 
The surface vector t is defined by 
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( 2 .«) 

For simp lie icy I it is assumed chat the fluid is a perfect gas wich 
constant molecular viscosity, thermal conductivity and specific heat. 
These assumptions are reasonable for low subsonic flows. By using these 
assumptions the right-hand side of equation (2.47) can be written as 
follows: 


7 • X ■ + 1/3UVB 


(2.49) 


Inserting equation (2.49) into equation (2.47) one gets 



(2.50) 


Also, 




(2.51) 


where n is the unit normal vector on the body surface, measured positive in 
the direction away from the solid surface. Because of the no-slip 

condition, V and — r— are zero everywhere on the surface. Thus, 

^ 3s^ 


7 




3u) 

3n 


(2.52) 
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Conbining the terns • equation (2*47) reduces to 

If . li |2 . (1/3)M II (2.53) 


If tiie pressure is non-^imensionalized with respect to the dynamic 
pressure at infinity, and all quantities are non-dimensionalized with 
respect to the reference quantities nentioned earlier, equation (2.33) 
becones 



2C 

Re '^3n 




(2.54) 


o 

where, C ■ — — and C is the chord length. 

2 * 

The dimensionless shear stress at the surface is given by; 


C 


f 



(2.55) 


knowing the surface pressure and the surface shear stress distributions, 
other quantities of interest such as lift, drag and moment can be easily 
obtained. 


Calculation of Loads 

Once the surface pressure and shear stress distributions are known, 
the loads are obtained from the following expressions. 
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S • ^ \ <2-56) 

C- - C- ♦ C_ (2.57) 

‘‘ *o *w 


C 


M 




where 




1 

C 


/ 


Cp‘9) % '*» 


(2.58) 


(2.59) 


(2.60) 


(2.61) 



(2.62) 


? I <=p<« {'ffi * ^ si} 
S'kc !’“<»> 


(2.63) 


(2.64) 
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The and C_ are force coefficients directed normal and tangential 
to the airfoil chord, and is the moment coefficient. The subscripts p 
and t denote the pressure and the skin friction contribution respectively. 
The moment is taken about the origin ot the coordinate system and is 
positive in the counterclockwise direction. 

The lift and drag coefficients referred to the wind axes are 
obtained from: 


• Cjj cos o- Qj sin a (2.65) 

Cjj ■ Cjj sina-«> cos a (2.66) 

where ^ is the angle of attack. 

Initial and Boundary Conditions 

The non-circulatory potential flow solution is used as an 
initial condition in Che present work. Along the body surface, the 
vanishing normal derivatives of enthalpy h and density p were used as 
boundary conditions for h and P. These conditions are convenient for an 
adiabatic wall. The surface values of the dilatation 0 were obtained using 
a three^point extrapolation formula during each iteration of the dila* 
tation transport equation. The boundary values of 8 vere relaxed and set 
to be zero vdienever the solution approached steady state. The integral 
expression (2.35) was used to compute the surface vorticity at each 
iteration while iterating the vorticity transport equation. 
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The potenciel flow veluet were weed ee in flow boundery condicionei 
«diile Che veniehing second deriveclves for h,P end cero vorcicicy were 
need es downecrena boundery condieione. The weke never epproeehed che 
downeereea boundery during che celculeCions. The e f or eaen Cloned boundery 
condieione preserve che ellipcic necure of the problem • 
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3. U80LT8 AKD DISCUSSION 


The procedure developed wee teeted on two probleae in order Co 
denonatrete Che ability of Che approach to coapute attached and aeparated 
flowa. The teat probleaa conaidered arat (i) laminar cn^reaaible flow 
around a circular cylinder i (ii) laadnar coapreaaible flow over an airfoil 
at aero angle of attack. Finally, the aethod waa applied Co the atatic 
atall problea. 

In the reaulta discuaaed below, the non-diaenaionaliaation ia done 
with reapecC to the free atreaa velocity and the characCeriatic length of 
the body in Che tranaforaed plana. In the airfoil caae and other teat 
caaea, the aolid body waa aet into motion iapulaively. Since the time 
rate of change of all flow variablea ia very high after the iapulaive 
atart, very aaall valuea of the time atep. At, are uaed to obtain proper 
Ciaewiae reaolution at the initial time levela. Aa the gradienta with 
reapect Co tiae decreaae, large valuea of AC are uaed. The under* 
relaxation paraaeter, which aoaetiaea controla the acceleration of the 
convergence of the iterationa, ia varied depending upon the type of problea 
conaidered. 


Vt^^ ^Bar_Caa preaaible Flow P eat a Circular Cylinder 

The preaent acheae alao haa been applied to the study of laminar 
coapreaaible flow paat a circular cylinder at a Reynolds niaber of 40, Mach 
number of 0.4 and Prandtl nuaber of 1. The Reynolds number is based on the 
cylinder diaaeter and the free stream velocity. This classical test case 
is chosen to demonstrate the ability of the approach to handle flows with 
massive separation. 


ORNUNAL PAGE IS 
OF POOR QUALITY 
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Tht grid aytcca oontiats of lia«a of conacanc radii and linaa of 
conacant angle 6. The linea of constant 8 are equally apaced with n /20 
intervals. In the radial direction, a stretching relation is adsuned as 
follows: 


r • e*, s • (j-l)As, j « l.Jaax 

By varying a uniformly, with ^ • 0.06, an exponential variation is 
obtained. The total number of grid points used is 2000 points. It should 
be noted that with Che present formulation not all of these grid points are 
involved in (he computations at all time levels. At the earlier time 
levels, the co^utacional region contained about 40Z of the total number of 
grid points. As the solution progressed in time, the number of points in 
the computational domain increased. When the computations were terminated 
at a time level of 15.1, che computational region contained all the 2000 
points. 

In order to co^ute the kinematic part of che problem che segmen- 
tation technique, explained in chapter II was used. The computational 
domain is divided into three annular regions R', R" and R"'. The inner 
region, R' , consists of 240 nodes. The intermediate region, R“, consists 
of 600 nodes, while the outer region, R'", contains 1160 nodal points. 
R^ions R', R” are matched a distance of .35 radii away from che 
surface; likewise the regions R" and R'" are matched at a distance of 2.525 
radii away from che eurface. The far-field boundary is located 17.916 
radii afay from the surface. The intagMl formuU la used Co 
compute che tangential velocity in region R* and on all of the boundaries. 
Then, the Poisson's equation is iterrted to get the velocities in R” and 
R'". 


ORIQINAL PAGE 18 
OF POOR QUALITY 
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Tht kiMtic «quation« (2.19)» (2.20), (2.21) «4i4 (2.22) «r« 

•pproxiaatcd hf an iaplicit fiBitc-dlffirMcc sehMa in th« poUr 
coordiiMtta and are aolved bj uaing the 'point aucceeeive under* 
relexetion* technique. Central differencea are uaed to epproxiaete the 
convection teraa. It ahould he noted that no ajaaetry wee eaauaed 
refarding thia prevent cave. 

The aolid body wee aet iuto aotion iapulaively. At thia iapulaive 
atert, the flow waa preacribed by the potentielflow aolution about e 
circular cylinder iaaeraed in a unifora atreaa. Tha tiae etep varied 

gradually froa 0*05 to 0*15. The aolution twa taminaced at t « 15.1. At 
thia tiae level the drag coefficient had converged to three digits. In the 
present caaet the tiae is non*diaena ions Used relative to the cylinder 
radius and the free atreaa velocity. 

In Figure 1 the surface pressure distribution at steady state is 
coapared with the ouaarieal aolution obtained by Sankar and Taaaa . 
The dgraeaant is quite good. 

In Table 1« the present coapreasible and incoapreaaible results are 
coeqtared with the cosiprcssible results of Reference 2 . In this Table, the 
separation angle d is aeasured froa the rear axis, and obtained as the 
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point on the surface where the vorticity changes sign. The length of the 
standing vortex (L/R) represents the distance between the center of the 
cylinder and the point on the centerline where the velocity changes sign. 
These conp^isons indicate chat the present solution and the finite* 
difference aethod give results that are in satisfactory agrceisenc. 


TABLE 1 


Comparison Between the Present Method and Reference 2 


• 

Presents 

Results 

Reference 2 


Incompressible 

Compressible 

Compressible 


IM=0) 

(M=0.4) 

(M»0.4) 

e 

sep 

53.5 

53.7 

52.5 

Lengch of Standing 

5.80 

5.85 

5.08 

Vortex (L/R) 




Pressure Drag 

1.025 

1.230 

1.32 

Coefficient Cj^ 




Friction Drag 

0.555 

.520 

.561 

Coefficient C. 




Total Drag 

1.580 

1.750 

1.881 

Coefficient C_ 

T 




Niniaua Surface 

Vorticity m . 

^ "min. 

-6.60 

-6.12 

-5.82 

Compressible 

Laminar Flow 

Past an Airfoil at 

Zero 


Angle of Attack 

The computational procedure developed here is next applied to the 
case of compressible laminar flow past a symmetric 9X chick Joukowski 
airfoil at a zero angle of attack* The airfoil is obtained by means of a 
conformal transformation of a unit circle* The chord leynolds number 
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considered in this case is 1000* The Mach nuaber is 0.4 and the Prandtl 
number is unity. The normalising reference time is obtained by dividing 
the transformed circular cylinder radius by the free steam velocity. All 
of the quantities are non-dimensionalised with respect to the free stream 
velocity and cylinder radius. 

The tangential velocity » , in the transformed plane, is 

calculated using the integral relation for the set of nodes on the first 
coordinate line next to the surface and at the outer boundaries. The 
Poisson's equation is then solved by using a 'successive point over* 
relaxation' technique in the rest of the domain. The difference kinetic 
equations, written in the transformed plane, are solved using a 'point 
under-relaxation' iterative technique. The circular domain is discretixed 
with 60 equally spaced points in the direction and 40 points in the 
radial direction. The time step is gradually varied from 0.0025 to 0.1 
during the course of the cosq>utations. The computations are initiated with 
an impulsive start. The initial surface vortex sheet strength is computed 
from the potential flow velocity values. The kineisatic computations are 
done with the finite Fourier series method. At a time level of 6.5, a 
steady state is determined to have been reached based upon the agreement 
(within IZ) of the computed surface vorticity values with those of the 
previous time level. 

In Figures 2, 3 the present surface pressure and surface vorticity 
values are cosq>ared with the corresponding values obtained in Reference 
2 . Both solutions are in very close agreement. The reference pressure 


used in these figures is the free stresa pressure. 


Lsainsr Coapressible Floe Psst an Airfoil st 
sn Angle of Attsck 

The sirfoil used in the present study is the 9Z thick synaetricsl 
JoukoHski airfoil described in appendix B. The chord Reynolds nuaber 
considered is 1000. The Mach nuaber is 0.4 1 the Prsndtl nuaber is 1.0 and 
the angle of attack is 15*^. 

A nuaber of publications (1, 2* 3, 5, 6, 7) have treated this 
problea before by incoii^ressible flow. It can be seen froa these results 
that the solutions are not quantitatively coaparable. However, there is a 
qualitative siailarity between the results. In the static stall case, the 
results depend on a nuaber of factors such as grid resolution, 
specification of the far-field boundary conditions and the nuaerical 
scheae. Therefore, in order to capture the saall coapressibility effects 
expected here, the coqtressible results have been oo^>ared with 

incoapressible results obtained using the saae coaputer prograa after 

'switching off the compressibility effects. 

Consequently, before solving Che coapressible static stall problea, 

it was appropriate to conduct first a series of incompressible numerical 
experiments Co : (i) test the code, (ii) inspect the sensitivity of this 

solution %rith the change of mesh size in the 6 direction, (iii) examine the 
role of Che ciae increment on the accuracy of the solution, (iv) examine 
the cyclic behavior of the solution, and (v) obtain incompressible data to 
be coapared later with the coapressible data. The difference between the 
t%ra solutions represents the effect of coapressibility. 
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Incowpwible Solution 

The iocoopressible solution for the ststic stall case was obtained 
by follAiii^ the saae procedure used later for the compressible case. In 
order to demonstrate the accuracy of the scheme, the incompressible 
solution has been cosqiared with the numerical results of Mehta C?]. As 
ahovm in Figure 4, the present results agree very well with Mehta's results 
at tike early time levels. As expected, the two solutions differ quanti- 
tatively at the later time levels. However, the qualitative behavior is 
similar at these later time levels. 

To illustrate the effect of the grid resolution on the solution, two 
sequences of solution were obtained for A9 » n/24 and A0 »7r/30. Figure 5 
shows the history of a load comparison between the two solutions. It is 
seen from this figure that the two solutions are comparable. Although 
there is no drastic difference between the two solutions, there still 
exists enough of a difference that there could be a misinterpretation of 
the results obtained for two different mesh size solutions, one compres- 
sible and die other incompressible. This experiment demonstrates the 
importance of using tiis same grid size whenever small cix^ressibility 
effects are examined. 

The continuation of the cyclic behavior of the solution for more 
than one cycle and the validity of the present isethod for a number of 
cycles of vortex shedding was demonstrated. The solution was advanced in 
time up to a dimensionless time level of 62 (the reference time being the 
transformed circle radius divided by the free-stream velocity). Figure 5 
shows the time history of loads which illustrate the cyclic b^avior of the 
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I 

[ 

|. solution with tine. Note that there are two cycles observed in the 

|* prescribed time range. This exercise provides considerable confidence in 

the fomulation of the problen and in the conputer prograa. 
j Finally, in order to study the effect of the tiae step on the 

solution, three nuaerical experiaents, with three different tiae 
I increaents, were perforaed. Table 2 shows the coq>arison8 aaong these 

I three solutions. Each solution has been started at tiae level of 20.175, 

and then advanced in tiae up to a tiae level of 21.615 and 24.735. The good 
I agreeaent aaong the solutions is apparent in Table 2. It could be 

concluded that, within the prescribed tine range, the sixe of the tiae step 

I ' 

plays a ninor role in the accuracy of the solution. 
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Table 2 


Coapariaon Aaong Three Different Time Increnents Solutions 


Tiae Level 


21.615 




S 


b) 

8, min 


0.2582 

0.1641 

0.1160 

-89.26 

0.2574 

0. 1654 

0.1150 

-89.41 

0.2623 

0. 1683 

0.1127 

-89.58 

0.3475 

0.1722 

0.1225 

-93.971 

0.3516 

0.1718 

0.1235 

-94.28 

0.3527 

0. 1686 

0.1257 

-94.46 


CoBpressible Solution 

lo the Figures that follow^ the chordwise distance denoted "chord 
percentage" is neasured from the leading edge of the airfoil. The force 
coefficients are normalized with respect to the free-stream velocity and 
the radius of the unit circle. The normal >.zed reference time is obtained 
by dividing the radius of the unit circle oy the free-^stream velocity. 

The grid system contains 48 equally spaced points in the 6-direction 
and 40 points in the radial direction. The exponential relation given in 
appendix E is applied for placing the points in the r direction. The time 
increment used in this numerical study is progressively increased from At * 
O.OOOS to At • 0.24. A total of 255 time steps were used to march the 
solution to a time level of 31.275, when the computations were terminated. 

Using the flowfield segmentation technique described earlier in 
chapter II, the velocity is obtained everyidiere in the computational 









I 
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doMln. Figure t shove the segmented compartments and the kinematic 
relation used in each of them. 

The iterative procedure used in solving the kinetic equations were 
varied to study their effects on the solution. Switching the iteration 
direction in the tangential coordinates was used to accelerate the con- 
vergence. The convergence criteria used in solving the vorticity transport 
equation was based on the maximum vortex strength variations between two 

consecutive iterations, fdiere the vortex sheet strength is defined by Y * 
2 

b)H dr. Invariably, the maximum variation between two consecutive iter- 
ations occured near the trailing edge, vdiich can be explained by examining 
equation (2.35). It is seen that the scale factor H appears in the 
denominator. Because the scale factor is very small near the trailing 
edge, it amplifies any error in the calculated value of the surface vortex 
sheet strength. The above criteria for convergence allows more tolerance 
for the vortex sheet strength near the trailing edge chan anywhere else. 
The tolerance level specified for the vortex sheet stren];.ch was 0.002 at 
Che earlier time levels and is subsequently reduced to 0.0005 at later time 
levels. Continuation of Che iteration beyond the above tolerance limit was 
not found effective in reducing Che residue. Hie residue instead oscil- 
lated around a minimum value without showing any tendency to reach zero. 

The maximum tolerance criterion used in iterating the enthalpy and 
density transport equations are 0.12 and .52 of the previous iteration, 
respectively. In iterating for the dilatation transport equation, a 
stringent tolerance limit in the vicinity of the airfoil was assumed in 
order to ensure proper convergence. In the outer regions and near the 
trailing edge, this limit is relaxed to accelerate the convergence. At 


later time levela, the oaxinun allowable tolerance, in the inner regionSi 
la taken to be 0.8Z. This represents the maxinua percentage variation 
between two consecutive iterations. 

As described earlier, the total number of nodes are 1920. However 
not all of these nodes were involved in the computations at all time 
levels. At early time levels, the vortical region is confined to only 

about 2SZ of the oaxiaium computational region. The computational time per 
time step depends on the extent of the computational boundary and varies 
from as little as 70 CPU seconds at the early time levels to 154 CPU seconds 
at later time levels on the CYBER'*70 computer with a C5C 6400 CPU. 

The average computational time required in the present study to 
advance the solution for one dimensionless time is 16 CPU minutes. Sankar 
and Tassa 2 used an ADI scheme to solve the primitive variable system of 
finite-difference equations, and took 11. S CPU minutes to advance the 
solution for one dimensionless time on the same computer. It should be 
emphasized here chat, in the present study, at later time levels Che memory 
requirements are larger than those required in Reference ^2. However, due 
to computer memory restrictions, the present computer program could not 
utilize Che maximum capacity of the CYBER-70 computer. Therefore, unneces- 
sary computations have been carried out for a number of time steps. 
Alternatively, if more computer core is used, along with using more 
sophisticated numerical procedures, it is believed chat the computational 
time requived to advance the solution for one dimensionless lime can be 
reduced below 11.5 CPU minutes. 

Table 3 gives the details of the time steps versus the computer time 


for the present computations 
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Table 3 


Sequence of Changing the Time Inctemenc and 
Sumary of Computer Tine 


No. of Time 
Steps 

At 

Time Level T 

Average CPU 
Time in Sec.* 

10 

0.0005 

0.005 

70 

10 

0.001 

0.015 

80 

10 

0.0015 

0.03 

85 

10 

0.0045 

0.075 

92 

10 

0.008 

0.155 

94 

10 

0.016 

0.315 

98 

10 

0.032 

0.635 

100 

10 

0.064 

1.275 

105 

20 

0.09 

3.075 

108 

20 

0.12 

5.475 

119 

30 

0.12 

9.075 

128 

25 

0.12 

12.075 

130 

80 

0.24 

31.275 

154 


* CYBER- 70/Model 74-6400 CPU. 

Flow Development 

The development of the flow field may be viewed as occur ing in four 
stages. These are: (i) impulsive start, (ii) formation and growth of the 
primary bubble, (iii) the bursting of the primary bubble which is 
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«$«oeiat«d with th« fonMtion of both the oecoadory «nd the trailing edge 
bubblee, end (iv) reettechaent of the priaery bubb7.e* The initial attached 
bubble which expands with tiae is referred as the 'priaary bubble*. 

The convention used is that the upper surface vorticity is negative 
for attached flows idiile positive vorticity indicates flow reversal. The 
opposite is true for the lower surface. 

The first stage of the flow field developaent reflects the effects 
of the iapulsive start. Inaediately after the iapulsive start, the 
vorticity is only non-zero at the surface, %ihile potential flow exists in 
the rest of the fluid. The rear stagnation point is located on the upper 
surface of the airfoil. Within a short tiae, the rear stagnation point 
moves close to the trailing edge. This aoveaent is associated with the 
foraation of a "starting vortex". At subsequent tiae levels the boundary 
layer starts growing on the upper and lower surfaces of the airfoil. The 
thickness of the boundary layer on both the upper and lower surfaces 
increases with tiae, as is observed froa the displaceaent of the 
'streaaline-like' lines near the surface. For convenience, the 
'streaaline like* lines will be called ' streaalines ' . The thickness of the 
boundary layer on the lower surface is saaller than the thickness of the 
boundary layer on the upper surface due to the existence of a favorable 
pressure gradient on aost of the lower surface. The extent of the region 
of adverse pressure gradient on the upper surface is shown in the preesure 
distribution plot. It is also observed, at this stage, that the aagnitude 
of the surface vorticity near the leading edge on the upper surface 
continues to decrease with tiae, forecasting the onset of separation in 
that nei^borhood when the surface vorticity changes sign. However, the 
separation does not actually occur until a tiae level of 1.88. During this 
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first st«fSy snd sftsr ths dsesy of the influence of the iapulsiee stertf 
the veloe of sterts inereeeing due to the growth of circulation after 
first reaching a ainiauai value at T ■ 1.0 as shown in Figure 7. The value 
of Gjj continually decreases because the decrease in the friction force on 
the lower surface as the positive vorticity decreases with tiaie. 

The second stage of flow developaent describes the occurrence of 
separationt together with the foraation and growth of the priaary bubble. 
The separation first takes place at 20X chord at T ■ 1.88. The sise of the 
separation bobble increases with tiae until it covers aost of the upper 
surface. This is expected* since the separation point aoves forward 
towards the leading edge and the reattachaent point aoves resrward towards 
the trailing edge. At tiae level 7.214* the separation snd ^e 
reattachaent points are about 95X chord length apart. Ae increases in the 
sise of the separation bubble increaaes the effective thickness of the 
airfoil* and the increase of the intensity of the reversed flow inside the 
bubble causes additional suction pressure on the upper surface. The above 
two factors result in an increase in the value of with tiae. During the 
duration of the priaary bubble* the drag coefficient reaains approxiaately 
constant. 

In the third stage of flow developaent* the priaary bubble is 
ruptured and an open bubble is foraed* indicating the cyclic atart of 
vortex shedding. The reattachaent point of the priaary bubble lifts off at 
a tiae level of 7.214 causing separated flow over alaost the entire upper 
surface. The increase in the nuaber of streaaline loops inside the 
separated bubble* indicates an intenaification of the reversed flow inside 
the bubble. The flow rotation inside the bubble is clockwise* with the 
fluid next to the surface aoving upstreaa. The preasure plots at 
T ■ 9.494* Figure 37* show a aaall region near the trailing edge idiere 
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th«re it a dacraaaa in the preasure in the direction of the aain flow 
outaide the bubble* Thia ia equivalent to an adverae preaaure gradient for 
Che flow near the aurface. A aaall counter clockwiae aeparation bubble 
appeara near the trailing edge at T ■ 11.775 aa a reault of the above 
mentioned preaaure gradient. The aize of thia bubble increaaea a lowly with 
time until it can be clearly aeen at T * 14.0SX. At thia time level, a 
aimilar adverae preaaure gradient developa at about S8X chordwiae diatance 
from the trailing edge. Thia reaulta in the appearance of a aecondary 
counterclockwiae bubble at T “ 16.71. The direction of the flow inaide 
thia bubble ia counterclockwiae, with the fluid near the aurface moving 
doimeCream toward the trailing edge. The intenaity of the flow rotation in 
the trailing edge bubble ia larger than it ia inaide the aecondary bubble, 
aa indicated by the number of atreamline loopa inaide chat bubble. The 
aize of the two bubblea increaaea with time. Aa time progreaaea, the 
primary bubble atarta to ahrink while Che other two email bubblea enlarge. 
The aecondary bubble expanda locally in the normal directi(m, whereaa the 
trailing edge bubble geta enlongated in the duwnetream direction. At thia 
atage of flow development, the lift coefficient keepa on increaaing, due to 
Che extent of the primary bubble beyond the trailing edge, until it reachea 
a maximum at T » 10.75. Heanidiile, the drag coefficient atarta increaaing 
very alowly after the time level 4.8. Thia alow increaae in the drag ia due 
Co Che increaae in the effective thickneaa of the airfoil aa judged by the 
chape of the zero atreamline. Thia cauaea an increaae in the preaaure 
drag. The value of continuoualy increaaea until it reachea a max i mum at 
T ■ 12. The downatream motion of the center of the ruptured clockwiae 
bubble, which ia accompanied by the appearance of the two amall counter** 
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clockwise bubbles, causes a t*o*'*i in the value of after 
T ■ 10.75. The reason for this drop is that the negative pressure 
sustained by the priaary bubble is partly reaoved by the foraiation of the 
two counterclockwise bubbles. ■ niis alao results in a decrease in the 
pressure drag, which explains the drop of after reaching a aMxianni at 
T - 12. 

The fourth stage of flow developaent involve the opening up of the 
secondary bubble, the lifting off of the trailing edge bubble, and the 
reattachaent of the upstream part of the primary bubble. The streamlines 
and equi-vorticity lines show the following flow development during this 
stage: (1) the secondary bubble splits the primary bubble and opens up to 
the outside flow at a time level of 20.51 { (2) the trailing edge bubble 
moves downstream until it lifts off the airfoil by T * 21.95; (3) the 
downstream part of the primary bubble starts to disengage from the surface 
at a tiS’^ level of 25.455; and (4) as time progresses, the upstream part of 
the primary bubble spreads in the downstream direction until the 
reattachment point reaches the trailing edge at a time level of 31.275, 
indicating the completion of the first cycle of vortex shedding. The 
streamline pattern at T ■ 31.275 looks similar to the pattern at the start 
of the cycle (T « 7.214), which indicates that the second cycle of 
oscillatory behavior is going to start at T * 31.275. As the secondary 
bubble opens up to the main stream, the reattachment point of the uputraam 
part of the primary bubble starts to move downstream, incraasing the region 
of the clockwise reversed flow. This reversed flow is able to sustain more 
suction pressure which results in an increase in the value of (^. The 
lift**off of the trailing edge bubble and the shedding of the downstream 
part of the primary bubble into the downstream flow enables the upstream 
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part of th« priaary bubblt to covar aoat of tha airfoil* Thia will load Co 
a continuoua iocraaaa in and aa Indicated in figuraa 8 and 9. 

Judging aithar froa tha contour plota of the atraaalinaa or froa tha 
tiea hiatoriaa of tha loada plota, it ia aatiaatad that ona cycle of vortex 
ahadding occura during tha tiaa period froa T ■ 7.214 to T ■ 31*279. With 
the airfoil chord aa the characteriatic length, the Strouhal niabar, 
defined C/(TV), idiara T ia tha period of the cycle, ia than found Co be 
0*195. 


Coapariaon Between tha Coapraaaible and Incoapreaeibla Solitiona 

In order to predict tha coapraaaibility affect for the praaant 
atatic atall caaa at a Mach niabar of 0.4, the coapraaaible and the 
Incoqtraaaibla aolutiona obtained uaing the aaaa grid are quantitatively 
compared in Table 4 * Tha iaportanca of uaing Che aaaM grid aixa in both 
aolutiona naa deaonatrated earlier in thia chapter* 

Baaed on the coapariaon ahown in Table 4, the obaervationa aade aay 
be auaaarixed aa followat 

(!) At the earlier tiae levela, the co^reaaibilicy aeeaa to 

decreaae Che rate of thickming of the boundary layer* 

(ii) The onaet of the aeparation of the priaary bubV>!e begina to 

appear at a later tiae level in the coapraaaible caae* 

(iii) Coapreaaibility aeeaa to play a aioor role in the growth of the 
priaary bubble* 

(iv) The co^reaaibility delaya the appearance of both the 

aecondary mid the trailing edge bubblea. Theae two bubblea 
grow at a faater rate in the coapreaaible caae than they grow 
in the incoapreaaible caae aa ahown in Figurea (8-10). 
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(v) When the firat cycle of vortex ehedding atarta, the difference 
in the force coefficienta between the coapreaaible and the 
inconpreaaible aolutiona ia snail. As tine advances, the 
difference between the two solutions gradually increases, in- 
dicating the increased influence of the conpressibility. At a 
tine level of approxinately 21 the conpressibility effect 
becones very snail. At this tine level, the tine rate of 
change of ilow variables decreases to a nininun. At later tine 
els, T 23 - 31.275, the coiapressibility effect appears to 
increase again but at a slower rate. 

(vi) The effect of conpressibility on the force coefficients is 
shown in Figures (8-10). This effect is conparable to the one 
computed in Reference 2, as shown in Figures 11 and 12. 
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Table 4 4 

Quantitative Comparison Between the Present Incompressible 
and Compressible Solutions 


Flow Events 

Compressible Data 

Incompressible Data 

Onset of separation of 
the primary bubble 

T - 1.88 

T » 1.76 

Separation location 
from leading edge 

20 Z 

23Z 

Cycle begins at T • 

7.214 

6.821 

Cycle ends at T * 

31.275 

29.815 

First appearance of 
trailing edge bubble 

T » 11.775 

T » 11.24 

First appearance of 
secondary bubble 

T » 16,70 

T = 15.374 

Opening up of the 
secondary bubble 

T « 20.51 

T = 19.81 

Strouhal number 
C/(TV ) 

CO 

0.15478 

0.16146 

max* 

1.32 

1.29 

niin. 

0.254 

0.262 


* The force coefficients are given after the recovery from the impulsive 
start. 


Comparison With Other Numerical Solutions 

The flow around a 9Z symmetric Joukowski airfoil at an angli> of 
attack of 15^, chord Reynolds number of 1000, Mach number of 0.4 and 

Prandtl number of one > has been solved numerically by Sankar and Taasa £23* 
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The conpuCational procedures of Reference 2 sre significantly 
different from those used in the present study, and it is very important to 
bear this in mind <dten comparing the results of the two studies. The 

procedures of Reference 2 are as follows: (1) the primative variables 

I 

(u,V,P,h) are used as the unknown flow variables; (2) the governing 
equations are discretized using central difference formulas for the 
spatial derivatives; (3) a second order artificial diffusion is added to 
the real diffusion term to stabilize the solution; (4) an AOl procedure ia 
used to solve the system of difference equations generated from the 
governing equations; (S) a fourth order dissipation term is added to the 
governing equations to eliminate the wiggles arising in the solutions; (6) 
the outer boundary is located at six chord lengths away from the airfoil 
(in Che present study the outer bounday is located at about 10 chord 
lengths away from the surface); and (7) a uniform flow is used to start the 
solution impulsively. 

Table 6 show typical comparison between the present results and the 
results obtained in Reference 15. 

The stability requirement for non-linear problems may impose more 
restrictions on the size of the time step even in the case of implicit 
schemes [s] However, Desideri et al [ 9 ] and Ballhaus ejt al* QoJ, in 
separate studies, have shown that a cyclic variation of the size of the 

f 

time step between two limits is helpful in obtaining convergence in the 
(AOI) schemes. In the present solution, no such restriction on the time 
step was required. It is believed that the use of the under*relaxation 
technique in solving the difference equations has a stabilizing effect on 
the solution. The maximum time step used in Reference 2 is 0.064, 
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whereas, in Che present study, the time size was successfully increased up 
Co 0.24. More details about the work can be found in Reference 11. 

Table 5 

Quantitative Comparison Between Che Present 
Results and Those of Reference 2 


Flow Feature 

Present Results 

Reference 2 

Onset of separation of 
Che primary bubble 

T ■ l.ft-^ 

T » 2.1719 

Separation location from 
leading edge 

2011 

24% 

Cycle begins at T * 

7.214 

6.3117 

Cycle ends at T ■ 

31.275 

27.376 

First appearance of Che 
trailing edge bubble 

T = 11.775 

T » 10.695 

First appearance of the 
secondary bubble 

T * 16.70 

T » 12.52 

Opening up of the 
secondary bubble 

T * 20.51 

T = 18.056 

Strouhal number C/(TV^} 

0.15478 

0.17626 


1.32 

1.34 

max. 



S . 

mxne 

0.254 

0.251 

s 

0.346 

0.364 

max. 



s . 

min* 

0.165 

0.141 
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4. THE USE OF THE INTEGRAL REPRESENTATION METHOD WITH 
SERIES SOLUTIONS FOR SOLVING THE NAVIER-STOKES EQUATIONS 


The use of orthogonal functions in solving the Navier-Stokes 
equations has offered high accuracy for certain probleas. The reason is the 
rapid decrease of the truncation error as the number of these functions 
used increase in a series representation of the solution (16). in this 
section of the report, the use of Fourier series with the integral 
representation method (1) is developed. The procedure is applied to simple 
test probleas. 

The Navier-Stokes equations, for steady incompressible flow, in a 
region R, with boundary B, can be written as follows (l): 


v(r) 




2v I .-,2 o 2g j 

-f, I 


(v X n )x(r - r) 
o o o 




dB 


2v J ,T :^,2 ■ 


dR 


,ss / ,, 

2 . B |. _.|2 o 

* o * 


-i- ^ ~ 


2ff B 




dB. 


dB. 


(4.1) 


(4.2) 


Where v and ware the velocity and the vorticity vector respectively, 'r is 
the posititm vector and n is the unit vector on B directed outwards. The 
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subscript o indicates that a variable or an integration is evaluated in the 
♦ 

space. Notice that (4.2) is nonlinear because it is equivalent to our 
familiar vorticity transport equation. 

In a polar (r,6 ) coordinate system, the vector equations (4.1) and 
(4.2) give the following scalar equations. 


, /Vo 

0 if : J ^ 


r sin(e - 0) 


r^ ♦ r^-2r r cos (0—6) 
o o o 


^v fr cos (6 -6)-rl 

^ r Lo o J 

^ R “5 7 

r r - 2r r cos(6 -0) 
o o o 


-v„ r sin(6 -6) 

1 f ®o ° ® 

■ ^ - 2r r oo.(S -0 ° 

o o o 


(4.3) 


, / OL r cos(6 -0) -r 

— ^ loo o 

^ J 


R r^ ♦ r^ - 2r r cos(0 -0) ® 

O GO 


♦ 4- ^ 


r 8in(0“0) 


r o 
o 


® + r^ - 2r r cos(0 -0) 

O GO 


, fej’^o 00.<e„-6) -r] 


- Jr r co«(6 -6) 

O 0 0 


(4.4) 


- r cos(6 -0)1 ♦ v„ <i) r sin(O-0 ) 
r : dS 


R - Jr r co.(G -0) 

O 0 0 


01 
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♦ _ J -2 _2 <IB 

® - 2r r cos (6 -0) 

o o o 


♦ 


1_ 

lit 


{%[’o - .. 

r-y y dB 

° - 2r r cos (0-0) 

o o o 


(4.5) 


In (4.3) - (4.5) and denote the velocity components in the r and 0 
direction respectively. Because we are looking at a periodic in the 0 - 
direction flowfield we can asatiae a solution in the form of finite Fourier 
series. 


N 

v^ » s^ + y (s^ cos n0 ♦ t sin n0) (4.6) 

r o n n 

n»l 


N 

2 (P- CO® “9 ”*■ ) (4.7) 

^ 0 _ , n n 

n»l 


N 

oi ■ a + i (a„ cos n0 sin n0) (4.8) 

o _• n n 

n»l 


In these equations, the Fourier coefficients ^n'*’ *o* 
s^'s, C^'®» Pg> dependent on r only. 

Using the method of residuals (Appendix c), the integrals in 
equations (4.3), (4.4) and (4.5) could be evaluated explicitly and only the 
Fourier coefficients are left to be determined (Appendix c)t 


s - 0 
o 


(4.9) 


, r r n+1 , 1 n-1 

*n • -5 / + ? J 

n 2 J n V oZ«inr c 


-2 •.<“ - 1 


(4.10) 
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- 2 J ■ ? J a (~) 

n z J n X o 2 J “o'r 


dr 


* i 'n<» » I P„(l) ' 


n-1 


r r 


o nL o r o 
o 

1 / 'o 1 A n-1 

"o ■ 2 / “„<r> “'o - 5 / «„<r> 


dr 


* i * 5 p„<i) 


n-1 


1 / 1 A n -1 

"n ■ 2 / 6„'r' “'o - 5 / e„<r> 


r o 


- 1 *.<>) '”■* ‘ 1 «„(») 


«« “ cu^l^ - Re f r dr 

O O *'^0 0 


o„ *a_(l) r 


n Re n /•■ , _ . 

-r' / <«„-«„) 


n n 


o J 
r dr 
o o 


* 2 / 

^ n n r o 

>« f* /r ,,r ,° 

• r J <«i. * "'o 

* o 

6„.8„(Ur"-f!r" / 


dr 


(4.11) 


(4.12) 


(4.13) 


(4.14) 


(4.15) 


(4.16) 


, Re r‘ x/^Ox“ 

r J \* <•*', 


47 


L ^ n n r o 
r 0 

where 1 <.n J^N, Che quantities knom 

velocity Fourier coefficients on boundary (or at r*l) and CtH > M sta 
Che coefficients from convective terms: 

N 

wv»C + 2 (L cos n6 + n sin **0 ) (4.18) 

r o u n 

n*l 

N 

u)v»M + y (w cos nS 4C sin n0 ) (4.19) 

o n *n 

n»l 

Once the Fourier coefficients of the velocities and the vorticity 
are determined, the velocity field and the vorticity field are easily cal- 
culated using equations (4.3) - (4.5). 

Equations (4.9)-(4.17) are solved using an iterative procedure for 
Che problem of steady flow inside a circle. The Fourier coefficients of 
Che boundary velocities are assumed to be known. 

Starting with known values of the Fourier coefficients of the 
vorticity , a^'s, superscript "i" being the iteration 

counter, the following steps constitute one iteration loop. 

(i) Determine the boundary values a^(l), a*^(l)'s and0*'(l)'s. 

on n 

The boundary values of Che velocities need to be satisfied by 

equations (4.9)-(4.14) during each iteration. Thus with one set of a^, 

a^'s and8*''s in R, the boundary values of a^(l), a^(l)'s and Sa'di's are 
an on n 

determined by 



o 


a r 
o o 


dr 

o 


p (1) 

*^o 


(4e20) 
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/ o <*r ■ P_U) * 

n o o *^ii n 


(4.21) 


/ «. C*‘ ’■'o ■ •»<" * ’.<*> 


(4.22) 


The equations (4.20)**(4.22) are the constrains on the boundary 

vorticity values derived from equations (4.9)-(4.14) when r«l. Tlie P^O)* 

p^(l)'s, q^(l)'s, s^(l)'s, t^(l)'s are the Fourier coefficients of the 

velocities on the boundary idiich are assumed to be known. With the proper 

numerical integration of the integrals in equations (4.20)-(4.22) , the 

values of a^'s and fi*''s at r>l are easily determined, 
o n n 

(ii) Compute the Fourier coefficients of the velocities in R. 

With the proper numerical integration, equations (4.9)**(4. 14) give 
explicit, point by point, calculation of ** ^o* ^n'** 

flow region R. 

(iii) Compute the Fourier coefficients of the convective terms. 

Because the Fourier coefficients of the velocities and of the 

vorticity are known at this stage, the Fourier coefficients C^,C*''s, 

on n 

li^'s, determined using equations (4.18) and (4.19). The 

coefficient need not to be determined because it will not get into the 
calculation of the Fourier coefficients of the vorticity. H.i8 quantity is 
associated with the static pressure level and it remains arbitrary when Che 
flow is incompressible. 

(iv) Cosq>ute the vorticity Fourier coefficients a a ^^^'s and 
8^*^'s in R. 

Equations (4.1S)-(4.17) permit explicit evaluation of a^^^'s 

o o 

«nd using quadratures if a^(l), a^(l)*s, 6^(1)'« and 

p^'s, r^'s are known on the rij^t hand side of these equations. 
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In the above iteration loop* it tns found necessary to employ a 
point under-relaxation technique to obtain converged solutions. The net. 
values of ® «nd 6 ^^^'s in R are computed from 


a ^ ^ ■ X a ♦ (1 -X )a*^ 
o o o o o 


(4.23) 


Xo ♦ (1 

n 


(4.24) 


X0 ♦ (1 

n 


(4.25) 


where 1 < n < N, X , X are the under-relaxation parameter and a*. a*'» and 
** *“ o on 

6^'s are the values computed in step (iv). 

Converged solution of a particular problem is assumed to exist idien 
the following criterion is satisfied. 


D • Max a . 
max ' oj 

J 


ft ‘ ft 1 a a ^ 

01., o. -o.,p. -p. 

OJ nj n j ’ nj n j 


<n < M 


(4.26) 


where subscript "j" denotes the Fourier coefficients at r«r.(r.< 1). A 

-4 

good value for e wss found to be 10 . 

The new numerical approach was tested and some of its features are 
discussed as the Reynolds number increases. 

(i) Plow without separation 

This is a closed streamline flow problem. The boundary velocities 


(4.27) 


so 




b 



6 


(4.28) 


Thia example %#aa treated by Burggraf (12). He uaed the Oaeen 
approximation, which took a rigid-body rotation as the basic flow. The 
solution of such linearised differential equation is only applicable to the 
case in which the whole flow field forms a singular circular eddy. The 
as 3 nsptotic vorticity value of the inviscid core at high Re he obtained is 
too low compared to Batchelor's (13) suggested model, which is proven to be 
quite adequate by Imai (14) and the present calculations. 

Using the present approach, converged solutions are obtained at 
different Reynolds numbers, from 0 to 1,000. The number of iterations, 
under-relaxation parameters and stream function values at the origin are 
presented it.’ Table 6 . The grid system is equivalent to 21 x 41 mesh 
points. The variation of maximum deviation D (in equation 4.22) versus 

IMX 

the iteration number is shotm in Figure 13 for Re ■ 1,000. In this case, the 
process of escalating the Reynolds number has not been used, i.e., the 
lower Reynolds number solution is not utilized as the initial solution to 
start the iteration procedure. This has otherwise been very effective in 
the calculations. 

The computer time for each iteration is about 1 second. Compared 
with the computer time that Imai had used at different Reynolds numbers 
(14} this approach showed speedup by a factor 2 to 3. 

The streamline patterns at Re ■ 0, SO, 300 are in Figure 14. The 
vorticity values at 0 ” 0 with different Re are presented in Figure IS. As 
Re increases, the value of vorticity in the inviscid core eventually goes 
to the value suggested by Batchelor and Wood (IS^l^. Also in this figure, 

* Wood actually calculated this value exactly after using the Von 

Mises transformation in the closed boundary layer region as Re'***‘” . 
Stated simply, the vorticity value in the inviscid core is 
determined by the root-sMmn-square speed of the closed boundary 
surface. 
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Ch€ boundary layar atruccura ia claarly aaan aa tha Rayoolda miabar 
incraaaaa. Figura 16 ahowa tha aitration of tha vortax cantar with 
incraaalng Rayitolda nuabar. It ia claarly aaan that tha ultiaate location 
of tha vortax cantar for Ra will ba tha cantar of tha circle. 

(ii) Plow conaiatina of two unavnatric racirculatinx raxiona 
Thia la alao a flow with cloaad atraaalinaa. Tha boundary 
conditiona are 


V • 0 (4.29) 

'b 

Vq ■ 3 ^ f (4.30) 

b 

Tha atreaalinea and vorticity valuaa at different Raynolda nuabara 
arc in Figure 17 and Figura 18. Notice that the aaynptotic flow pattern for 
hi^ Raynolda nuabera cannot be found in a aiaple way aa in the laat 
exaaplei aince the fora of the aeparation line ia net known in advance. 

(iii) Onflow - outflow problea 

The flow problea together with the boundary velocitiea are depicted 
in Figure 19 In thia caae, the flow ia ayaMtric about the x-axia» thua the 
Fourier expreaaion could be aiaplifiad to 

N 

hi ■ y 6 ain nO (4.31) 

n 


V ■ 

r 


N > 

y a coa n6 
n-1 " 


(4.32) 


S2 


«nd 



«v. 


wv^ - 


9 

y o aio nd 
n-1 ^ 

R 

X sin n6 
n-1 “ 

N 

n * X cog n0 

I n 


(4.33) 

(4.34) 

(4.35) 


Due to Che velocity dlecontinuity et the boundary, e larger number 
of grid pointa ia uaed in the agrimuthal direction. The grid ayaten uaed 
in thia caae ia equivalent to 21 x 81 meah pointa in the whole plane. The 
81 puinta along the circuaiferential direction were proven adequate. The 
uae of 121 pointa in that direction generated differencea of the 
atreamf unction valuea leaa than 1 percent compared with Che reaulta that 
uaed 81 pointa. The effect of the grid aiae in radial direction haa alao 
been teated, aa ahown in Table 7 . The computer time for each iteration, on 
the 21 X 81 grid ayatem, ia about 1.3 aeconda, a CDC 6600. 

The atreanline patterna with different Reynolda numbera are 
depicted in Pijture 20. The occurrence and growth of the aeparation bubble 
are .ilearly aeen in thia figure, aa the Reynolda number increaaea. 



Table 7 


Re 

Ar 


ITER 




O 


^HAX 

100 

.1 

.2 

47 

.111 


.05 

.1 

80 

.123 

200 

.1 

.06 

130 

.116 


.05 

.05 

211 

.132 

300 

.1 

.03 

230 

.115 


.05 

.002 

440 

.131 
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RGURE(3 ) JOUKOWSKI AIRFOIL AT ZERO ANGLE OF ATTAG 


SURFACE VORTICITY OISTRIBUTION AT STEADY STATE 









Region I The Integral Relation (3.30) is used to Calculate 

Region II The Integral Relation (3.30) is used to Calculate 

Region III FinitcDifference Equation (3.34) is used to Calculate V 


FIGURE ( 6) SEGMENTATION OF KINEMATICS AT LATE TIME LEVELS 
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FIGURE( 8 ) TIME HISTORIES Of LOADS 




FIGUREC 9 ). 




FIGURE(IO) TIME HISTORIES OF LOADS 




FIGURE(ii) time histories OF LOADS 













Figure 16 
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tfPENDIX A 

OBUVATXOR OP THE RXinTXC TRANSPORT EQUATIONS 

(1) Tb< Vorticity Trampoft Equation 

Hm vortieiCy tr«n«port •quaeion (2.17) say b« wrieten «• 


♦ ij <A.l) 


For CMO-diamtional flowt wt have 


- 1 ; <f K <m|* • M |j) I7 <»■ I7- ju«) 


8 /I 8 /«„ 8u 2 , a\ ‘ 8 , 3u . 8v. 


(A.2) 


Vx(wu) ■ - (V.7w«u6) • -7.(^w) 


(A.3) 


Upon diffaraneiacing tba tana la aquation (A>2) it reducaa to 

7.T) " Tj+Ij (A.4) 

tdiara T| rapraaanea cha auaution of tha fraa-viacosity changa taraa and T 2 
capraaanta Cha Catma which ineluda tha change of viacoait^. 
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H If - l!i;> ‘ If 




♦ M 




^ ? “*y 


• fl ^^*“y *\^x' “x * f <*y®^ 


♦ M (- I A ♦ I 

y P ^ 


fl ^yrv 


- -4 (2P V -u p„ - y. 


For siaplicicy and conveniancci the viacoaiCy coefficienta |i and 
X(X ■ -ju), apeeific heat radio Y and the theraal conductivity k have bean 
eonaidarad to be eonatant in the praaant ai^aonic atudy. Xn ganeralt the 
incorporation of variable coafficianta haa no conceptual affect on the 
•athod. Several nunarical axpariaanta have bean parfomed to verify the 
foragoiof aaauaption and the raaulta indicate a negligible contribution of 
the non-conatant coefficient tana 

By uaing the above aaau^itionaf equation (A«6) reduces to the 
following aiaple font 


aut 

n 


- ♦ ( 




9(^*ByP«h^ 


U.7) 


idieret 






♦ j • — ~ (y (fcx 7P).k - 


. Y-l ♦ ♦ ♦ 

^2* y * ^**'*‘‘ 


i« Cht unit vxetor notMl eo Uit plant of flow (x*y) 


(2) Ibo Diltttcioo Trtntoort Ituttioo 


Iho dllaettion crtntporc oquacion (2-18) at/ bt wriecan at 


• - 5.(V.^)V - 2^ 7^ - 2^ h 7^1oP - 2^ ^loP.7h 


1 ^ 1 ♦ ♦ 

♦ fe V4 7.T) 


U.8) 


Sptcialiting chit for Cvo-diatnaional flow and txpanding Cht firtc 
and cht late cars in cht righc-hand tidtf ont obcaint 


•v.(v.5)» . - |j (u |S . V |a> - 2y (u U . » |2) 


* *" a? * ^ ■ ^ * *^* * ’ Isy’ 


(A.9) 


Furchart nocing chac* 
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B > iii 4 ^ 

ax ay 


Equation (A*9)taduca« to Che following fora 


♦ ♦ 36 as 2 2 

“^(V.V)V ■ -(u ■^♦V' 5 - + u ♦v +2v u) 
ox oy X y X y 


Equation (A. 10) can be further reduced to 


-7.(V.7)y “ ♦ 8 ♦ 27vx7u.k) 


- - v7.(V© ♦ 27vx7u.k) 


Siailarly, 


•}; (f < 1 ; (a* I; -f uB) *57 (u 57 *m |;»> 




By separating Che right-hand side into two terms, similar 
has been done in the vorticiCy transport equation, one ubtains 


■* 1 

7 .(i 7 .t) • C, ♦ C- 
P 12 


(A. 10) 

(A. 11) 

(A. 12) 
to tdiat 

(A. 13) 


where 
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A u 2 A U ♦ U ± 

Cj • I J 7^B - (j ^ 7P.7B ♦ ^ (7P*7w).k 


a. 14) 


Cj - (| P) 76.^ ♦ I (W5p).k - (|p)87^P 


♦ (4) (u +vU ) + (v+u)m 

p XXX yyy x yxy 


- (^)(V, * Uy) ♦ Py .p,) . (|) 1 5 p.5« 




(A. 15) 


Ncglecciog Che cent C^t for che reasons mentioned before in (I)* Che 
dilacacion Cransporc equation is written as follows 


II » - 7.(VB) ♦ (37^) 7.^6 ♦ X<“fB»P»h) 

where 

X - Xi ♦ X2 ♦ X3 
Xj^ ■ ”2 (^v X 7u).k 

X, - — |- 47 p *?0 ♦ (^px^).kj 

p Re ' •* 

X3 • - (^) (V^h ♦ hV^lnp ♦ ^.^np) 


(A. 16) 






I r 

[i 


•3 


(3) Th» ItwriT Trwoort Equation 

Otiat Ch« •quation of state (2.4) in the energy equation (2.3), one 
obtains the following 

Iy ■ -V.^ - (Y-l)6h - J 7.q ♦ J (VV:T) (A.17) 

Equation (A. 17) can be written in diaensionless fora as follows 

H • - "-“B" - 5 ‘ sn; <*•»*> 

Expanding the conductive and dissapative teras in equation (A*1E) 
one gets the following expressions 


PRe.Pr 


Y 1 3 


3h 


V.a • i — s — ^ — (-k i::-) 
P Re.Pr 3Xj ^ j 517^ 


U.19) 


and 


Y 1 Y ^**1 


(A.20) 


Specializing these two equations for two-diaensional flow, they 
reduce to: 




PRe.Pr 


(A.21) 


and 




M«gl«ctitig tiie themal conductivity variation tana and placing 


aquation (A.21) (A*22) into aquation (A*18) one gets; 

|~ “ ♦ 9(w,p,B,h) 


^ara 


0j - (2-K)0h 

®2 * ♦ 4(Vv X Vu).k| 


(4) Tha Density Equation 

Tha continuity aquation (2.2) can be written as follows. 


1 dP 

f at 


V 

♦ ^ . 7p 
P 


Further 9 noting that 


0 « 7.V 


(A.23) 


1A.24) 


thus aquation (A. 24) baconea, 
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• -6 -V.^lnP 


(A.2S) 


la coaaarvativa fonif aquatioa (A<25) can be wrieten as followa 


1^ - -5.(VlnP) ♦ r 


(A.26) 


eherc 


r - 3(lnP-l) 
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APPEMDIX 8 

T86IISF0RM6T10N RBUTX0N8 

1h« transfoxMtion relationa uMd in cht derivation of the kinetic 
tranaport equationa era given here. The following expressions relate tiia 
aMeheaatical operatior's done in the phjrsical planS) ph, with thoaa that 
perforaad in the tranaforaed plane* T* 


(Vf.V 

ph. 


(7fxV^) 

ph. 


(T^f) 

ph. 


(7. (fV)) 


IT T 


- ~ (7f X %) 

T 


- ^ (T^f ) 
H T 


• ~ (7. (f^)) 
T 


(B.l) 


(B.2) 


(B.3) 


(B.4)' 


where* 

(i) R is the transf onset ion scale factor 
(ii) f and g are scalar functions which are invariant with the 


transformation 
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(iii) it tilt vtlocity vector in the phyticnl plane, while it 

the apparent velocity vector in the trantfortied plane and ita 
eoaponantt (V«,V ) are given by equationa (2*27) and (2-28). 

tf V 

The coaponentt of the velocity vector V in the phytical plane (u,v) 

« 

are related to the velocity coaponenta follova: 

u a ~ (-V ^ — V (B.S) 

" g2 ' ^8 dr r ’'r d8' 


V 



r dr 


♦ 


Ip iL) 

r de' 


(B.6) 


where y aid x are the Carteaian coordinatea in the phytical plane. The 
traaaforaad quantitiea can be expanded in polar coordinatea (t-8) at 
follaaa: 


H • If ♦ 7 H • M] 


(B.7) 


‘ <H • If - H • lf> 


(B.8) 


2 .2. 




i£ ♦ i iii (B.9) 


T 3a 


2 'r(r-c) ^r-c^ ^ ^ r 


where a ia the stretched radial coordinate in cb** transformed polar plme 
and ia given by; 


8 ^ 

r * e ♦ c 


(B.IO) 
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APP88PIX C 

EVALOATIOR OF INTECRAL8 IN BQDATIOHS (4.3) - (4.5) 
The iotegrals that need co be evaluated are 


J - r 
m o 


r 

m o J 


^ 2 

J - r' 
m o 


/ 


2 2 , 
r * T 

= a d 

Ir r cos(0 -0) 

0 

o o 

sin m6 

r cos (6 “0)-r 

o 

^ o J 

r^+r^-2r r cos (6 -0) 

o 

0 o 

cos m6 

sin(0 -0) 

0 

0 

r^+r^-2r r cos (8 -0) 
o o o 

sin m 6 

sin(0 -0) 

o 

® d 


(C.l) 


dd 


dO 


r^+r^-2r r coa(d -S) 
o o o 


(C.2) 


(C.3) 


(C.4) 


where 0^ N. If Che following complex variablea are defined as 


W ■ I + i J 

mm m 


(C.5) 


W • I + i J 
mm m 


(C.6) 


and 


ie 

z a e ■ cos 6 + i sin 6 

o o o 


(C.7) 


i6 « • . « 

a ■ e ■ COS0 1 sin6 


(C.8) 
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iihere i The quanticiea and could be reexpreaead as follows 

after the variable transformation from 6 » 8 to St s t 


W . 1 f!2L__l!o 

V (* - I_ 


♦ zh- 2rs. 
o 




dz^ 


or o 

o 


a) 


m-1 f 2 2. 

(s -a ) 
o 


W - ! 

m 2t Y 

J (*_ - I- 


da 

o r 'X'o-r** 
o 


(C.9) 


(C.IO) 


By using of residue theorem, Che integrals in equations (C.9> and (C.IO) 
are evaluated explicitly and the results are 


(1) for r < r : 
o 


W ■ 0 (C.ll) 

0 

m-1 

W a® (C.12) 

“ ''o 

• 0 (C.13) 

m-1 

W ■ (“) ** (C.14) 

ID r 

o 

(2) for r> r : 
o 

r 

W ■ - 2ir(— ) (C.15) 

o r 

r^ m+1 

K- ■ -»f<— ) a <C. 16 ) 

Is r 

W - 0 (C.17) 

O 

^ r m+1 

a® (C.18) 

ID r 
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I 


<t 

r 

> 

y 

I 

5 . 

) 


[ 

H 

r 

1 ; 

i: 

1 

1 

I 

! 

1 


mhers 1 < m < N. The values of I t J 1 Z 1 

* **“ B n B 

by equations (C.5) end (C.6). Thus 

J^'s could be determined easily 

(1) for r < r : 
— 0 


0 

1 

0 

(C.19) 

^ m-1 


I ■ iK— “/ cos m6 

« '0 

(C.20) 

j ■ 0 
0 

(C.21) 



r 

J ■ ir(--) sin m6 

® '0 

(C.22) 

9H 

0 

1 

0 

(C.23) 

m-1 


■ - ir(^) sin me 

“ >^o 

(C.24) 

^o"0 

(C.25) 

B*1 

^ r 

J • iK-“) cos me 

(C.26) 


O 


I 

1 


i: 

1 ; 

[ 

[ 



(2) for r> r : 
o 

r 

Iq • - 2lK^) (C.27) 

r «+l 

I - -ir (7^) CCS me (C.28) 

ID r 

Jj, - 0 (C.29) 

r Btfl 

J ■ -ir (~) sin ffl (C.30) 

01 r 

Ig • 0 (C.31) 

r sH>l 

■ - ir(-^) sin me (C.32) 

01 r 

“ 0 (C.33) 

r m+1 

J_*ir( 7^) cos m0 (C.34) 

01 r 
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The integrele thet need to be evelueted in kinetiee ere 


I- • » 

ai o 


o 4 - 2r r cos(9 -0) 

o o o 


de 


(C.35) 


J - r 


/ 


irein m0 , r - r coi(6 - ® | 
o L o ' 0 J 


d6 


o - 2r r cos(6 -0) ° 

o 0 o 


(C.36) 


I - r r / 
m o J 


2ir coe o0 ein(0 -0) 

O 0 


d0. 


o r ♦ r “ 2r r coe(0 -60 

O 0 0 


(C.37) 


J - r r / 
■ 0 J 


2ir sin sin(0 -0) 
0 o 


<*e_ 


o r + r -2r r cos(0 -0) 

O O 0 


(C.38) 


also the coaplex quantities are 


W • 1 4 i J 


(C.39) 


W - I 4 i J 
am a 


(C.40) 


Following the same procedure as in kinematics, the following expressions 
ere obtained 


a 2ir 




r-: • 




<*. - f - r •> *'* 

o 


(C.41) 




[ 

K 

[ 

r 

r. 

i; 

i; 


•-1 /_2 . 2 , 


5 .i/!2_i!sl£L_ 




<•. - r •>“. - r •> 

o 


and the final reaulta are 


(1) for r < r I 
o 

1 ■ 2ir 

o 


■ w(^)" coa b6 

o 

•'o"® 

"W (““)* *in ■ 
o 

T ■ - ir(--)" ain a6 

■ ir(~)" ain ad 
IB r 

o 


(2) for r > r^j 
o 


I • 0 
o 

r m 

I ■ - ir(r^) coa a6 

u V 

J " 0 
o 

r a 

J ■ -ir<^) ain ad 




I ■ -ir (~) Bin me 
B r 

T • 0 

o 

^ r a 

■ ir(— ) coa a6 
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(C.43) 

(C.44) 

(C.45) 

(C.46) 

(C.47) 

(C.48) 

(C.49) 

(C.50) 


(C.Sl) 

(C.52) 

(C.53) 

(C.54) 

(C.55) 

(C.56) 

(C.57) 

(C.5«) 


where 1 < a < N 
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itfPBIIDIXD 

OBRZVATIOII OF MLATtOlfSIlIPS BBTtfBBN FOORXBR OOBmCXOITS 


The following oxprnooion for io obtained after aubatituting 
aquationa(4.6, 4.8 ) into equation (4.3) 

’r - f. /[«. K ‘ j, <Vn 

♦ 2 I- ♦ J«> 

ZirpLo o n n n nj 

n*i 

- JiU*" K ‘ ", <«•.<» *, * V*’ 3.>] (».» 

n*i 

where It It Jt J ere thoae quantitiea defined in Section (i)t Appendix B. 
Froa equations (C.ll) > (C.?4) 


1 J B r r n+1 r n*l 1 

■ » { •*" “* ■ *”* **J 


dr 


1 1 nr r r 1 

- -r j i 1®. ainn0- 6(“) coa ooldr 

^ V n-lL ” 'o * 'o J * 


2 n.iL“ » 

♦ ^ Z |pg,<l) »i“ n0 - q^(l) r***^ coa nsj (0.2) 


If is also exp^dad in Fourier aariaa (aa in aquation 4.4) )t the 
corresponding Fourier coefficients are related by 


♦ I .„(!) r“-‘ 


- i q„(l) 


(16 


where 1 < n 


( 4 iich ie 


1 ^ fk ®^1 . I n«f 

'» • - 5 / “.<r> - i /. .(J-) 

* J ‘„( 1 > r"-* . ^ ,.-l 


dr 


( 17 ) 


< H. SiaiUrly, the expreeeion for froa equation (4.4) 

''• * 5 / [», 'o * <«.I„ ♦ e„ J„)] dr. 

* I„ ♦ t^(l) T^)] 


n»l 

M 

2 

n«l 

N 

I 

n»l 


^ 2ir ^0 ^ 5 <Pn(l) !„ ♦ q„(l) J )] (D. 3 ) 


1 / ?r r n*l , 

on-1- ^ «‘>« ne ♦ B„(^) aionejdr^ 


’e-i/ sCvr 

o n»l 

1 P r a*l Q j 

■ ' i “• "® ‘ -J-ir 

O 

* 5 "« - ‘r<» «« .«] 

* ? 2 «•• «9 ♦ « (I) r"'* .1. .al 


/n A \ 
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Expanding in Fourier aeries as in equation ( 4,7) the correapondiog 


coefficients are related by equations (4.12)“(4.14). 


By expanding a,Vg in equations (4.ig) (4.19) equation (4.5) could 


be expressed as 


w “ - J CC I ♦ 1 (£ I ♦ n J )] dr 
2 v J ‘•^o o n n ’n n •* < 

o n*i 




n*i 


(D. 5 ) 


where I, It Jt J's are those quantities defined in Section (ii) Appendix 


Ct and f t f * g 'a are the Fourier coefficients of the total pressure, 
o* n’ *n 


h « f + 2 (f cos n6 + g sin n$ ) 
° n-1 " " 


(D.6) 


By the results of equations (c.43)-(c>58)t equation (D.5) becoaes 


w ■ ^ / 2 (~) “® * ^^n * nsjdr^ 


o n»l 


-r /' K* Vf->° 

♦ - Mjj) sin nojjdr^ 


♦ I* I [- yu t- .in .e . yi) r“ «>. «•] 


6-Z 


i 
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R 

♦ * ♦ i o (1) r" cos 00 + B (1) sin n6 (D.7) 

n-1 " “ 


Thus I in using the Fourier series expansion for b) (equation 4.8 ), the 
corresponding coefficients are related by 


a ■ J ‘**“0 

00 y ® ^ 


/ 


( 21 ) 


»» ■ r / «n - “'o - r / <S. ♦ «.xf >" "'c 




(D.8) 


»« • ?/ ‘'V *“ ><r)‘ * r / <•‘0 - ".xf >“ ^ 


f (1) r" ♦ 7B„(1) r" 
2 n 2 n 


(D.9) 


idiere 1 _< n ^ R. Equations (d> 8) and (b*9) are equivalent to equations 
(4.16) and(4.i7) respectively, where the terms g^^O) <nd f^d) are cancelled 
by applying the equations (d* 8) and (D.9) on the boundary. 
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APPBIIDIX B 

THB INTBGRAL RBPRB8BRTAT10N8 FOR A DOUBLY-COMIIBCTBD RB6I0II 

The boundary integrals in equations (4.3)and (4.4)conaist of two 
parts* and B 2 




(i.i) 


1 “2 


At the interior boundary B 2 * the velocity coaponents are assuaed to 
be known. The Fourier coefficients of the velocities at B 2 are aero except 
p^(R>* which is the aagnitude of the circuaferential velocity at r^ due to 
the solid rotation. 

The equations VD.l) for the velocity v^ needs to be rewritten to 
include the B 2 effect: 




* h 'n * '.<*> •’,>3 


- 7 - Cs„(R) 1^*1 («„(R) ♦ t„(R) J )] 

ZIT o o n n n n 


- k k<‘> ‘ 


n 

* k h * J, ‘ (i.2) 


where t^CR)* 9^(9) 1 Fourier coefficients 






f 

6 
r 
r 
r 
r 

f r 
r 
r 
r 
r 
r 
r 

i: 

r 

r 

[ 
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of th« velocitUi at r-«. ihu. 


•o " 0 


(15) 


1 i* “♦i J 

"if dt^ * j 8n<r>“’^ 


* i - 1. _n-l 

« n 


dr 


!%<•> 


(S.3> 


1 o+l 1 „_, 

‘"■■5^“n<r> "'o-i /««<r> 


dr 


1 

2 "o' 


‘ i t.(l) r"-‘ . 1 , (,) 
A n 


(E.4) 


where 1 < n < N. Similarly, the velocity ^ 


i./i. 


18 reexpressed as: 


* S [•„<» ?o • j, <•.(«) \ * t^d) ?.)] 
■ f. «o * (•„(«) I„ ♦ t^(«) J__)] 


n«l 

- f.tp.<«) ?. * (P,<*) i„ . ,^<«) T_,] 


(B.S) 


obtained 


»y -.in* equtiM. (c.l7) Cl» follo»ii„ 


relations are 



irtier* *nd g^CK) tht Fourier eooffieimto 

prosouro «t r«B. By uoing the equetione (C.43)-<C.58)t 


.r R r n 


« ■ I* jJ [<«„ - Cb>“» ♦ (1, ♦ M 


- r / CW“* 

♦ 2 [- f|,<i)r” sin n6 ♦ g^(l) r“ co# nO] 


n»l 


R 


- 2 [- C(»)<7)" •in ne ♦ g^CRXf) coe nO] 

2 n-1 ” ' ^ ' 


of the totel 

^)elo u0]dr^ 
n o 

)ein nOldr^ 

i O 


(B.IO) 
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